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Fast  roll  control  in  fruit  flies:  evidence  for  a  nonlinear  response  mechanism 


Tsevi  Beatus1  John  M.  Guckenheimer2,  and  Itai  Cohen1 
Departments  of  1  Physics  and  2 Mathematics,  Cornell  University,  Ithaca,  New  York  14853,  USA 

(Dated:  November  10,  2014) 

Statement  of  the  problem  studied: 

Due  to  aerodynamic  instabilities,  stable  flapping  flight  requires  ever-present  fast  corrective  actions. 
Here  we  investigate  how  flies  control  body  roll  angle,  which  is  unstable  and  their  most  sensitive 
degree  of  freedom.  We  glue  a  magnet  to  each  fly  and  apply  a  short  magnetic  pulse  that  rolls  it 
in  mid-air.  Fast  video  shows  flies  correct  perturbations  up  to  100°  within  30  ±  7ms  by  applying  a 
stroke-amplitude  asymmetry  that  is  well  described  by  a  linear  PI  controller.  For  more  aggressive 
perturbations,  we  show  evidence  for  nonlinear  and  hierarchical  control  mechanisms.  Flies  respond  to 
roll  perturbations  within  5ms,  making  this  correction  reflex  one  of  the  fastest  in  the  animal  kingdom. 


Summary  of  the  most  important  results: 

Locomoting  organisms  evolved  mechanisms  to  control 
their  motion  and  maintain  stability  against  mechanical 
disturbances.  The  control  challenge  is  prominent  in  small 
flying  insects  since  their  small  moment  of  inertia  renders 
them  susceptible  even  to  gentle  air  currents  [1—4] .  More¬ 
over,  they  fly  at  Reynolds  numbers  Re  =  102  —  104,  in 
which  flows  are  unsteady  [5,  6].  Most  importantly,  flap¬ 
ping  flight  is  aerodynanrically  unstable,  on  a  time  scale 
of  a  few  wing-beats  [7  18].  It  is,  therefore,  intriguing 
how  insects  overcome  such  control  challenges  and  man¬ 
age  to  fly  with  impressive  stability,  maneuverability  and 
robustness,  outmaneuvering  any  man-made  flying  device. 

Among  the  body  Euler  angles  -  yaw,  pitch  and  roll 
roll  is  most  sensitive  to  perturbing  torques  since  the  mo¬ 
ment  of  inertia  of  the  insect’s  long  axis  is  smallest  [1,  2], 
Recent  fluid  dynamics  simulations  suggest  roll  is  unsta¬ 
ble  due  to  an  unsteady  aerodynamic  mechanism,  where  it 
is  positively  coupled  to  sideways  motion  via  asymmetry 
of  the  leading-edge  vortex  attached  to  each  wing  [14-18]. 
Such  results  indicate  flies  can  lose  their  body  attitude  due 
to  roll  perturbations  within  4  wing-beats.  Controlling 
roll  is  also  crucial  for  maintaining  direction  and  altitude. 
Thus,  any  basic  understanding  of  insect  flight  demands 
quantitative  analysis  of  roll  control. 

Previous  studies  used  tethered  animals  to  measure  re¬ 
sponse  to  imposed  roll  rotations  [19-24]  and  visual  roll 
stimuli  [23-29].  In  such  experiments,  however,  the  teth¬ 
ered  insect  does  not  control  its  motion  and  often  exhibits 
wing  kinematics  and  torques  qualitatively  different  from 
those  in  free-flight  [30,  31].  More  recently,  free-flight  ex¬ 
periments  used  vortices  [1-3]  and  impulsive  gusts  [4]  to 
perturb  insects,  highlighting  the  sensitivity  of  the  roll  an¬ 
gle  to  perturbations.  Understanding  roll  control  however, 
requires  fast  and  accurate  quantitative  measurements  of 
wing  and  body  kinematics  in  response  to  controlled  mid¬ 
air  perturbation  impulses  -  a  methodology  recently  ap¬ 
plied  to  study  yaw  control  [32] .  Crucially,  these  previous 
works  typically  consider  only  the  linear  response  [29,  32- 
38].  Whether  non-linear  mechanisms  come  into  play  in 
natural  free  flight,  where  both  large  and  coupled  pertur¬ 
bations  are  common  [1]  remains  unknown. 


Here,  we  perturb  a  fruit-fly  ( Drosophila  melanogaster ) 
by  gluing  a  magnet  to  its  back  and  applying  a  ~  5ms 
magnetic  pulse  that  rolls  it  in  mid-air.  The  perturbation 
amplitude  ranges  from  roll  deflections  of  5°  to  multiple 
rotations  along  roll.  We  use  high  speed  video  to  film  the 
fly’s  corrective  maneuver  and  measure  its  wing  and  body 
kinematics  [39] .  We  find  that  for  roll  perturbations  up  to 
70°  flies  generate  corrective  torques  by  applying  a  stroke- 
amplitude  asymmetry  that  is  described  by  the  output 
of  a  linear  PI  controller  for  roll.  The  asymmetry  starts 
only  one  wing-beat  (~  5ms)  after  the  perturbation  onset, 
making  the  roll  correction  reflex  one  of  the  fastest  in  the 
animal  kingdom  [40].  Surprisingly,  however,  we  find  that 
linear  control  is  not  sufficient  to  explain  the  response  for 
multiple  rotations  nor  the  overarching  control  structure 
for  simultaneously  handling  yaw,  pitch,  and  roll. 

To  exert  mid-air  mechanical  perturbations  along  roll 
we  glue  a  magnetic  pin,  1.5  —  2mm  long,  to  the  dor¬ 
sal  thoracic  surface  of  each  fly  (Fig.  1).  The  pin  is  ori¬ 
ented  horizontally  and  perpendicular  to  the  body  axis. 
In  each  experiment  ~  15  flies  are  released  in  a  transpar¬ 
ent  chamber  equipped  with  two  Helmholtz  coils  that  are 
used  to  generate  a  vertical  magnetic  field  (~  10_2Tesla). 
When  a  fly  crosses  the  filming  volume,  a  laser-trigger 
initiates  video  recording  at  8000  frames  s^1  along  three 
orthogonal  axes,  as  well  as  a  5ms  (1  wing  beat)  mag¬ 
netic  pulse  [32,  39].  Since  fruit  flies  fly  with  their  body 
axis  pitched  up  at  ~  45°  and  since  the  moment  of  iner¬ 
tia  along  their  body  axis  is  ~  4  times  smaller  than  the 
other  axes,  the  largest  deflection  is  generated  along  the 
body  roll  axis,  with  smaller  perturbations  along  pitch 
and  yaw  (Fig.  la,b).  Using  a  custom  image  analysis  al¬ 
gorithm  [39],  we  extract  a  3D  kinematic  description  of 
the  fly  (Fig.  1)  consisting  of  its  body  position  and  ori¬ 
entation  (Fig.  lb)  as  well  as  the  Euler  angles  (Fig.  lc) 
for  both  wings  (see  Supplementary  Information  (SI)).  We 
analyzed  20  sequences  that  span  a  perturbation  range  be¬ 
tween  5°  and  100°,  in  which  the  flies  perform  a  steady 
flight  before  and  after  the  correction  maneuver. 

A  representative  example  of  a  fly  recovering  from  a 
60°  roll  perturbation  is  shown  in  Fig.  1  and  Movie  1. 
Body  Euler  angles,  roll  (p),  yaw  (</>b)  and  pitch  ( 6 b)  are 
plotted  in  Fig.  le.  The  magnetic  field  was  applied  be- 
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FIG.  1.  Roll  perturbation  and  correction:  (a)  Images  from 
3  orthogonal  cameras  of  a  roll  perturbation  and  correction 
maneuver.  Each  panel  shows  4  superposed  images.  The  3D- 
rendered  fly  represents  the  kinematic  data  of  the  body  and 
wings.  The  perturbation  location  (red  line)  is  shown  on  the 
fly’s  center-of-mass  trajectory  (green).  In  the  second  snapshot 
the  fly  is  rolled  60°  to  its  left,  (b)  Definition  of  body  Euler 
angles  with  respect  to  the  lab  frame.  Xb  is  the  long  body  axis, 
(c)  Definition  of  Euler  angles  and  body  frame  (xb,yb,Zb). 
Wing  angles  are  measured  in  the  body  frame  with  respect  to 
the  stroke  plane  (shaded  blue,  SI),  (d)  Top  and  side-views  of 
10  consecutive  wing-strokes  of  the  maneuver,  taken  when  the 
wings  are  at  their  forward-most  position.  The  perturbation 
wing-beat  is  numbered  0.  (e)  The  body  Euler  angles  during 
the  maneuver.  The  perturbation  was  applied  between  0  —  5ms 
(yellow).  White  and  gray  stripes  represent  forward  and  back 
strokes,  respectively.  Yaw  and  pitch  were  sampled  at  8000Hz. 
Roll  was  measured  manually  at  the  middle  of  each  half-stroke 
and  smoothed  by  a  spline  (dashed  line).  Measurement  errors 
are  comparable  to  the  symbols  size. 


tween  t  =  0  —  5ms  and  induced  a  maximum  roll  veloc¬ 
ity  of  7000°s_1  resulting  in  a  deflection  of  60°  within 
t  =  13.5ms  (Frame  3,  Fig.  Id).  The  fly  recovered  its  ini¬ 
tial  roll  angle  within  35ms,  or  8  wing-beats.  Top  views 
show  a  clear  asymmetry  in  wing  stroke  angles  during  the 
maneuver,  which  starts  a  single  wing-beat  (5ms)  after 
the  perturbation  (Fig.  Id,  frames  1  —  4).  During  the 
maneuver  the  left  wing  stroke  amplitude  increases  while 


the  right  wing  stroke  amplitude  decreases.  The  fly  also 
spreads  its  legs  from  their  folded  position  (frames  4—8)  as 
in  a  typical  landing  response  [41-43].  In  addition,  smaller 
deflections  of  25°  left  in  yaw  and  5°  down  in  pitch  were  in¬ 
duced,  since  the  applied  torque  is  not  completely  aligned 
with  a  principal  body  axis  (Fig.  le,  SI). 

We  quantify  the  asymmetry  in  wing  kinematics  by 
plotting  the  wing  stroke  angles  during  the  maneuver 
(Fig.  2a).  We  find  large  differences  -  up  to  70°  -  between 
their  peak-to-peak  amplitudes  (Fig.  2b,  SI).  The  ampli¬ 
tude  asymmetry  began  1  wing-beat  after  the  onset  of  the 
perturbation,  and  lasted  for  5  wing-beats.  The  flapping 
frequency  of  both  wings  remained  nearly  constant  during 
the  maneuver.  Hence,  to  maintain  the  amplitude  asym¬ 
metry,  the  right  wing  moved  faster  than  the  left  (Fig.  2c). 
To  first  order,  this  difference  in  velocity  leads  to  asym¬ 
metry  in  the  aerodynamic  forces  of  the  two  wings  and 
generates  a  correcting  torque. 

To  calculate  the  aerodynamic  torque  generated  by  the 
insect,  we  used  the  full  measured  wing  and  body  kine¬ 
matics  combined  with  a  quasi-steady-state  model  for  the 
aerodynamic  force  [44]  produced  by  each  wing  (SI).  The 
calculated  torques  are  similar  for  other  quasi-steady-state 
force  models  [45,  46]  as  well.  The  components  of  the  aero¬ 
dynamic  torque  vector  along  the  Xb  and  Zb  body  axes 
were  averaged  over  half-strokes  and  plotted  in  Fig.  2d,e 
(see  Fig.  lc  for  axes  definition).  The  torque  magnitude 
is  roughly  5nN  •  m  and  is  comparable  to  torques  exerted 
by  tethered  fruit  flies  [23].  Both  the  Xb  and  Zb  torque 
components  have  a  corrective  effect  along  roll,  and  both 
exhibit  distinct  peaks  (solid  circles)  that  appear  simulta¬ 
neously  with  the  stroke-amplitude  asymmetry.  Our  cal¬ 
culations  also  show  that  passive  roll  damping  is  too  weak 
to  explain  the  observed  maneuver  (SI). 

This  amplitude  asymmetry  can  be  described  by  the 
response  of  a  linear,  proportional-integral  (PI)  controller: 


A<F(f)  =  Kpp(t  -  AT)  +  Kip(t  -  AT). 


(1) 


Here,  the  output  A<f>  is  the  difference  between  the  right 
and  left  wing  stroke  amplitudes,  and  the  controller’s  in¬ 
put  is  the  body  roll  velocity,  p,  which  flies  measure  using 
their  gyroscopic  sensor  system  associated  with  the  hal- 
tere  organs  [20,  47,  48].  The  controller  is  defined  by  three 
parameters:  the  proportional  gain  Kp ,  the  integral  gain 
K[,  and  a  delay  AT  that  describes  the  neuro-muscular 
response  time.  Fitting  these  three  parameters  using  the 
measured  p,  p  and  A$  (Fig.  2f),  we  find  this  controller 
response  (red  curve)  is  sufficient  to  reproduce  the  time- 
dependence  of  A<f>  (green  circles).  Moreover,  the  fast 
rise  time  can  be  attributed  to  the  term  proportional  to 
the  roll  velocity  (blue  curve).  Simpler  models,  such  as 
I-  and  P-controllers,  can  be  excluded  (SI).  As  such,  this 
PI  model  is  the  simplest  continuous  linear  control  mech¬ 
anism  consistent  with  our  observations. 

The  salient  correction  feature  -  a  wing  stroke  ampli¬ 
tude  asymmetry  with  a  delay  of  ~  1  wing-beat  -  was 
observed  in  all  the  recorded  events  (Fig.  3a).  Moreover, 
this  feature  was  robust  to  variability  in  initial  flight  pose 
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FIG.  2.  Roll  correction  mechanism  for  the  correction  ma¬ 
neuver  shown  in  Fig.  1.  (a-c)  Wing  stroke  kinematics  versus 
time:  (a)  The  stroke  angle  tj>  of  the  right  (red)  and  left  (blue) 
wings;  (b)  their  peak-to-peak  amplitude  <F,  and  (c)  their  an¬ 
gular  velocity  (j>  (c).  (d-e)  Mean  aerodynamic  torque  along 
each  half  stroke,  calculated  from  the  measured  wing  kinemat¬ 
ics  using  a  quasi-steady  state  aerodynamic  force  model.  Solid 
symbols  highlight  the  correcting  wing  strokes,  (d)  the  torque 
component  along  the  body  axis  Xb-  Positive  torque  induces 
a  corrective  right  roll;  (e)  The  torque  component  along  Zb. 
Negative  values  have  corrective  effect,  (f)  Wing  stroke  am¬ 
plitude  difference  A<F  (green),  and  a  fit  for  a  PI  controller 
(eq.  1,  red),  with  AT  =  4.4  ±  0.25ms,  Kp  =  6  ±  0.5ms,  and 
Ki  =  0.7  ±  0.05  (fitted  values  ±  confidence  intervals,  shaded 
red,  SI).  Contributions  of  the  1st  and  2nd  terms  of  Eq.  1  are 
shown  in  blue  and  black,  respectively.  Measurement  errors  in 
(a-f)  are  comparable  to  the  symbols  size. 


Ap  [deg]  Max  roll  velocity  [deg  s'1] 


FIG.  3.  (a)  Response  time  histograms  for  multiple  perturba¬ 
tion  events:  Delay  time  of  fitted  PI  controllers  (dark  gray) 
and  the  time  to  reach  A<f>  =  10°  measured  from  wing  kine¬ 
matics  (green).  Histograms  are  shifted  to  improve  visibility, 
(b)  Mean  roll  acceleration  versus  A<f>  for  299  forward  strokes 
(red)  and  304  back  strokes  (green).  The  correlation  coeffi¬ 
cient  of  these  two  quantities  is  —0.67  (p-value  <  1CF12).  The 
dashed  line  has  a  slope  of  —20  ■  103s-2.  (c)  Roll  correction 
time  Tc  measured  in  wing-beats  as  a  function  of  A p  for  mul¬ 
tiple  events.  The  mean  correction  time  (thick  dashed  line)  is 
6.8±1.6  wing-beats  (meanistandard  deviation).  Thin  dashed 
lines  indicate  la  margins,  (d)  Maximum  roll  displacement  Ap 
as  a  function  of  the  maximum  roll  angular  velocity  for  multi¬ 
ple  events.  The  dashed  line  has  a  fitted  slope  of  9.9ms  with  a 
95%  confidence  interval  of  ±0.75ms. 


and  velocity  (Movies  2,3)  and  observed  even  for  two  con¬ 
secutive  perturbation  pulses  (Movie  4,  Figs.  S2,3).  In- 
depth  analysis  of  11  correction  maneuvers  showed  the 
Ad>  response  is  consistent  with  the  PI  controller  model 
(fit  mean  R 2  =  0.86).  Fitting  control  parameters  for 
each  maneuver  separately,  we  find  Kp  =  4.8  ±  2.4ms  and 
K[  =  0.6  ±  0.3  (meanistandard  deviation).  The  mean 
response  time  AT  =  4.6  ±  1ms  is  comparable  to  a  single 
wing-beat  period  (Fig.  3a).  The  confidence  intervals  for 
parameters  in  each  event  are  smaller  than  standard  de¬ 
viation,  indicating  differences  between  fitted  controllers 
can  be  attributed  to  natural  variation  between  flies  (SI). 

This  response  time  is  3.5  times  faster  than  the  flies’ 
response  to  yaw  perturbations  [32]  and  2.5  times  faster 
than  their  response  to  pitch  perturbations  [13].  The  roll 
response  is  also  5  times  faster  than  visual  startle  response 
in  flies  (25ms)  [49]  and  is  completed  before  they  can  elicit 
a  visual  response.  In  fact,  this  ability  to  respond  within 
5ms,  places  this  reflex  among  the  fastest  in  the  animal 
kingdom,  comparable  to  those  in  blowflies  (5  —  7.5ms) 
[24],  teleost  fish  (5  —  10ms)  [40,  50],  trap-jaw  ants  (4  — 
10ms)  [51]  and  cockroaches  (11  —  18ms)  [52,  53]. 

To  further  illustrate  that  wing  asymmetry  generates 
corrective  roll  dynamics,  we  determine  the  mean  roll  ac¬ 
celeration  for  ~  600  half-strokes  and  plot  it  versus  Ad) 
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(Fig.  3b,  SI).  The  A4>  asymmetry  is  negatively  correlated 
with  the  roll  acceleration  generated  by  the  fly.  Thus, 
for  example,  negative  A<!>  is  correlated  with  positive  roll 
acceleration  (Figs.  1,2,S2,S3).  The  data  also  show  two 
hallmarks  of  linear  control.  First,  the  correction  time  is 
insensitive  to  the  maximum  roll  deflection,  A p.  Here,  the 
correction  time  Tc  is  the  time  between  onset  of  the  pertur¬ 
bation  and  the  moment  when  the  roll  angle  reaches  10% 
of  A p.  Plotting  Tc  as  a  function  of  A p  shows  the  correc¬ 
tion  time  is  6.8±1.6  wing  beats  (mean±SD,  n  =  20)  with 
little  dependence  on  perturbation  amplitude  (Fig.  3c). 
Second,  we  find  A p  increases  linearly  with  maximum  roll 
velocity  (Fig.  3d).  Collectively,  these  data  suggest  that 
as  with  yaw  [32],  the  response  to  roll  perturbations  is 
well  described  by  a  reduced  order  model  of  a  linear  PI 
controller  with  time  delay. 

To  test  the  linear  control  model  we  challenged  the  flies 
with  extreme  perturbations  in  which  they  were  spun  mul¬ 
tiple  times  in  mid  air  by  a  series  of  magnetic  pulses. 
The  fly  shown  in  Fig.  4a  and  Movie  5  was  rotated  8 
times  to  its  right.  The  accumulated  roll  angle  exceeded 
3000°  (Fig.  4b)  and  the  maximum  roll  velocity  was  over 
60,000°s_1  (Fig.  4c).  During  the  perturbation,  the  fly 
was  unable  to  oppose  the  magnetic  torque.  We  cap¬ 
tured  3  such  events,  all  showing  the  same  behavior.  Re¬ 
markably,  once  the  magnetic  pulses  stopped,  the  flies  re¬ 
gained  control  within  3  —  4  wing-beats.  Our  calculations 
show  roll  deceleration  is  only  explained  by  active  flap¬ 
ping,  rather  than  passive  damping  due  to  the  wings  (SI) . 

If  the  fly’s  roll  controller  were  a  linear  PI  (Eq.  1),  the 
integral  term  should  have  accumulated  a  signal  corre¬ 
sponding  to  a  3000°  deflection.  The  resulting  correction 
maneuver  would  require  the  fly  to  rotate  3000°  in  the 
opposite  direction.  Clearly  flies  circumvent  this  scenario 
by  employing  something  like  an  integral  “anti-windup” 
operation  that  prevents  such  accumulations  from  taking 
place.  It  is  plausible  that  in  such  maneuvers  the  time 
scales  are  long  enough  that  flies  incorporate  an  additional 
sensory  modality  to  determine  the  direction  of  gravity 
and  do  not  rely  on  integration  of  angular  velocity.  Either 
way,  the  observed  behavior  is  an  example  of  a  nonlinear 
feature  for  roll  control  during  extreme  perturbations. 

We  observed  an  additional  nonlinear  control  feature 
that  arose  because  the  perturbing  torque  was  not  per¬ 
fectly  aligned  with  a  principle  body  axis  and  induced  sec¬ 
ondary  deflections  along  the  yaw  and  body  velocity  direc¬ 
tions.  In  contrast  to  previous  experiments  where  pertur¬ 
bations  were  restricted  to  yaw  and  were  corrected  with 
small  errors  [32],  here  we  observed  that  the  secondary 
deflections  along  the  yaw  and  body  velocity  orientations 
were  often  left  partially  uncorrected  (Figs.  4de).  For  ex¬ 
ample,  even  after  the  roll  correction,  deflections  at  t  =  11 
wingbeats  show  yaw  “errors”  larger  than  10%  of  the  ini¬ 
tial  yaw  perturbation  (above  the  green  region  in  (Fig. 
4d)  and  body  velocity  orientation  “errors”  that  are  even 
larger  than  the  initial  body  velocity  orientation  perturba¬ 
tion  (orange  area  in  Figs.  4e).  The  qualitatively  different 
yaw  response  with  and  without  the  roll  deflection  hints 
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FIG.  4.  Extreme  perturbation,  (a)  Snapshots,  taken  1.25ms 
apart,  of  a  fly  being  rotated  at  ~  60,  000°s_1  by  a  pulse  train 
of  magnetic  torques,  (b)  Roll  kinematics  during  the  maneu¬ 
ver  (blue  circles),  and  its  smoothing  spline  (black  line).  The 
perturbation  was  on  between  t  =  0  —  72ms.  Measurement  er¬ 
rors  are  smaller  than  the  symbol  size,  (c)  Roll  velocity  versus 
time  (black  line).  The  fly  was  unable  to  control  its  flight  dur¬ 
ing  the  perturbation  (orange  shading)  and  actively  corrected 
roll  once  the  perturbation  stopped  (green  shading).  The  fly 
leveled  itself  to  p  =  20°  within  4  wing-beats  after  the  pertur¬ 
bation  stopped,  towards  p  =  0.  (d)  Yaw  difference  between 
t  =  0  and  t  =  11  wing-beats  as  a  function  of  the  maximum 
yaw  perturbation,  (e)  A  similar  plot  for  the  velocity  heading. 


that  the  fly  prioritizes  roll  correction  while  compromising 
other  angles.  Such  prioritization  would  require  a  control 
mechanism  with  nonlinear  coupling  of  the  body  angles. 

Collectively,  these  results  shed  light  on  important  as¬ 
pects  of  the  fly  controller.  First,  the  flies’  response 
to  extreme  perturbations  of  multiple  roll  rotations  in¬ 
dicates  they  employ  a  nonlinear  mechanism  such  as  in¬ 
tegral  “anti-windup”.  Second,  their  response  to  com¬ 
bined  perturbations  along  multiple  axes  suggests  their 
coupled  control  mechanism  is  hierarchical  and  therefore 
nonlinear,  prioritizing  roll  -  an  unstable  degree  of  free¬ 
dom  -  while  compromising  other  angles.  Consequently, 
mapping  out  the  architecture  and  hierarchy  of  the  full 
controller  will  necessitate  measurement  of  the  insect’s 
response  to  sophisticated  perturbations  along  multiple 
axes,  whose  direction,  amplitude  and  timing  are  individ¬ 
ually  controlled.  Finally,  the  rapid  response  time  of  5ms 
underscores  the  importance  of  roll  control  for  flight.  In 
addition,  this  small  time  scale  suggests  the  neural  cir¬ 
cuit  governing  the  response  must  be  simple,  consisting 
of  only  a  few  neural  synapses  between  the  haltere  organs 
and  the  steering  muscles.  We  envision  that  this  simplic¬ 
ity  opens  the  door  to  combined  bottom  up  neuroscience 
experiments  and  top  down  biolocomotion  studies,  such 
as  those  presented  here,  aimed  at  elucidating  the  func- 
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tion  of  each  neuron  in  the  control  circuit  of  this  important 
model  animal.  Such  measurements  will  further  reveal  the 
neural  architecture  behind  the  elegant  strategies  insects 
use  to  manage  their  actuation  resources  and  achieve  the 
grace  and  performance  of  their  flight. 
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